Fermi surface of the filled-skutterudite superconductor LaRu 4 Pi 2 : A clue to the origin 

of the metal-insulator transition in PrRu 4 P 12 
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We report the de Haas-van Alphen (dHvA) effect and magneto-resistance in the filled-skutterudite 
superconductor LaRu4Pi2, which is a reference material of PrRu4Pi2 that exhibits a metal-insulator 
(M-I) transition at Tmi — 60 K. The observed dHvA branches for the main Fermi surface (FS) are 
well explained by the band-structure calculation, using the full potential linearized augmented-plane- 
wave method with the local-density approximation, suggesting a nesting instability with q = (1, 0, 0) 
in the main multiply connected FS as expected also in PrRu4Pi2. Observed cyclotron effective 
masses of (2.6 — 11.8)rao, which are roughly twice the calculated masses, indicate the large mass 
enhancement even in the La-skutterudites. Comparing the FS between LaRu4Pi2 and PrRu4Pi2, an 
essential role of c — / hybridization cooperating with the FS nesting in driving the the M-I transition 
in PrRu4Pi2 has been clarified. 

PACS numbers: 75.20.Hr, 75.30.Mb, 71.18.+y, 71.27. +a 



The filled-skutterudite compounds i?Tr 4 Pni2 (R = 
rare earth, Tr= Fe, Ru, Os; and Pn=pnictogen) have 
attracted much attention^ for exhibiting novel phys- 
ical properties, i.e., metal-insulator (M-I) transition, 2 
unconventional superconductivity^^ semiconducting be- 
havior.""' magnetic ordering^ and their prospect in ther- 
moelectric applications. 1 Particularly interesting are the 
Pr-based skutterudites in which quadrupolar interactions 
are believed to play an important role for the anomalous 
propertiesi^SiL&S PrRu4Pi2 shows an exotic M-I tran- 
sition at Tmi — 60 probably with a non-magnetic 
origin, as suggested by no anomaly in the magnetic sus- 
ceptibility (x) and no essential effect of the magnetic field 
on the specific-heat anomaly up to 12 Although 
several studies have been carried ou t)2iiS4iiii24iii4 the 
mechanism of the M-I transition remains one of the most 
mysterious puzzles in the skutterudites. 

From the band-structure calculation(i^i2iil Harima 
has suggested that the M-I transition could be triggered 
by the nearly perfect nesting of the main Fermi surface 
(FS) with q= (1,0,0); the shape of the main FS is nearly 
cubic and its volume is one-half of the body-centered cu- 
bic (bcc) Brillouin zone (BZ). Actually, both the elec- 
tron 11 and the x-ray 12-13 diffraction experiments have re- 
vealed a doubling of the unit-cell volume across Tmi; it is 
suggested that the structure changes from bcc (Im3) to 
simple cubic (Pm3), caused by the different shift of Ru 
and P ions around two adjacent Pr ions, i.e, between the 
corner and body-centered Pr ion sites in bcc structure. 
However, the predicted distortion of P-ion to obtain an 
insulating band structure is much larger than that ex- 
perimentally determined by the x-ray diffraction^SiiLiS 
If the 4/ electrons in PrRu4Pi2 are localized, as indi- 
cated by the Curie- Weiss behavior of xr 1 or only weakly 
hybridize with conduction electrons, then the electronic 
structure is expected to be basically the same as the ref- 
erence material LaRu4Pi2, which shows a superconduct- 



ing transition at Tq ~ 7 Kji£ without a 4/ electron. In 
that case, LaRu4Pi2 is also expected to have the same 
nesting instability, however no M-I transition has been 
reported. Thus, it is essentially important to clarify the 
FS of both compounds experimentally to understand the 
origin of the M-I transition in PrRu4Pi2. However, the 
high transition temperature of Tmi — 60 K prevents the 
observation of dHvA oscillations in the metallic state of 
PrRu4Pi2. As an alternative way, the experimental de- 
termination of FS in LaRii4Pi2 should be useful to un- 
derstand the origin of the M-I transition in PrRu4Pi2, 
if it is carefully compared with the computed FSs of the 
two compounds based on the band-structure calculation. 

We have succeeded in growing high-quality single crys- 
tals of LaRu4Pi2 and in observing the dHvA effect, the 
preliminary result of which has been briefly reported in 
Ref. 20. In this paper, we report the detailed experi- 
mental results on the dHvA effect and magnetoresistance 
along with the band-structure calculation in LaRu4Pi2, 
followed by a discussion on the possible origin of the M-I 
transition in PrRu4Pi2. 

Single crystals of LaRu4Pi2 were grown by the tin-flux 
method which is basically the same as described in Ref. 
5. The raw materials were 4N (99.99% pure) -La, -Ru, 
6N-P, and 5N-Sn. The residual resistivity ratio (RRR) 
of the present sample is about 700, indicating the high 
sample quality. The dHvA experiments were performed 
in a top loading 3 He cryostat down to 0.4 K with a 15 T 
superconducting (SC) magnet and a top loading dilu- 
tion refrigerator cooled down to ~ 30 mK with a 17 T 
SC magnet. The dHvA signals were detected by means 
of the conventional field modulation method with a low 
frequency (/ ~ 10 Hz). The magnetoresistance was mea- 
sured by the dc four-probe method using a top loading 
3 He cryostat equipped with a 16 T SC magnet. 

Figurenshows (a) a typical recorder trace of the dHvA 
oscillations and (b) the corresponding fast Fourier trans- 



formation (FFT) spectra in LaRu4Pi2 at 31 mK for the 
magnetic field (H) along the (111) direction. The dHvA 
signals have been observed just above the superconduct- 
ing upper critical field Hq2 ~ 36.5 kOe, indicating also 
the high quality of the sample. Figurc^a) shows the an- 
gular dependences of dHvA frequencies both in the (010) 
plane and the (110) plane. The calculated dHvA frequen- 
cies are shown in Fig. ffib) , the description of which is 
given later. In the experiment, six fundamental dHvA 
branches labeled a, /3, u>, S, e, and p have been ob- 
served. The largest frequency branch a is observed over 
the whole field angles both in the (010) plane and (110) 
plane, although it is not perfectly continuous. The j3, 
u), and p branches are observed in the limited angular 
range around [111] direction. The 5 branch is observed 
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FIG. 1: (a) Typical dHvA oscillations and (b) its fast Fourier 
transformation (FFT) spectrum in LaRu4Pi2. 
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FIG. 2: Angular dependence of the (a) experimental and 
(b) theoretical dHvA frequencies in LaRu4Pi2. The dHvA 
branches indicated by the closed marks in the Fig. EJb) are 
identified by the experiments. 



(48th-band) hole FS in LaRu4Pi2, which is similar to 
that of LaFe4Pi2- LaRu4Pi2 is an uncompensated metal 
with a different carrier concentration of electrons and 
holes, since it has one molecule per unit cell. If TMR 
increases with the field, then an open orbit exists for the 
direction perpendicular to both the current and field. In 
the present experiment for H\\ [100], TMR increases with 
H 11 , suggesting the existence of an open orbit. 
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in a narrow angular range centered at the [001] direction. 
The e branch is observed in a wider angular range in the 
(010) plane except near [100]. 

We have measured the transverse magnetoresistance 
(TMR), Ap/p = [p(H)-p(0)]/p(0), in LaRu 4 Pi 2 to con- 
firm the existence of open orbit on the FS. Figure[3]shows 
the field dependence of TMR for the field along the [100] 
and [101] directions and the inset shows the angular de- 
pendences of TMR under the constant fields in the (010) 
plane. As was also observed in LaFe4Pi2r^ the angular 
dependence of TMR in LaRu4Pi2 is highly anisotropic, 
suggesting the existence of an open orbit. Actually, as 
shown later, we have observed the multiply connected 



FIG. 3: Magnetic field dependence of the transverse magne- 
toresistance (TMR) for H along the [100] and [101] directions. 
The inset shows the angular dependence of TMR under the 
constant fields of 7 T and 14 T in the (010) plane. 



In order to assign the origin of dHvA branches, the 
band-structure calculation has been performed by us- 
ing the full potential linearized augmented-plane-wave 
(FLAPW) method with the local-density approximation 
(LDA). We used the room-temperature lattice constant 
a = 8.0608A and the fractional coordinates of P at the 
24g site as x = 0, y = 0.3594, and z = 0.1434 for the cal- 
culation^ The details of the calculation are essentially 
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the same as in Ref. [H Figures 0] and [3] show the cal- 
culated energy band structure and the FS in LaRu4Pi2, 
respectively. The FS is composed of the 46th, 47th and 
48th band hole sheets. The main difference with the FS 
of PrRu4Pi2 is the existence of the 46th and 47th small 
hole sheets in LaRu4Pi2, while PrRu4Pi2 has only a mul- 
tiply connected 49th band hole-FS sheet with one band 
counted for two localized 4/ electrons^ The 46th and 
47th bands form nearly spherical sheets centered at the 
T point, which stretch slightly along the (111) and (100) 
directions, respectively. The 48th band gives a multi- 
ply connected hole-FS sheet centered at the T point. Its 
shape is roughly a round cube whose volume is nearly 
one half of the BZ size and has a good nesting vector q= 
(1,0,0) as expected also in the 49th-band FS of PrRu4Pi2. 
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FIG. 4: Energy band structure of LaRu4Pi2. 

The calculated dHvA frequencies, shown in Fig. OJb), 
reasonably reproduce the angular dependence of the ob- 
served dHvA branches a, (3, 5, e, and /x except uj. The a' 
branch is found to originate from the electronlike orbits 
in the 48th band FS. The other dHvA branches predicted 
in the band-structure calculation have not been observed 
in the present experiments, and the probable origin of 
this discrepancy is discussed later. We have also esti- 
mated the cyclotron effective mass m* from the temper- 
ature dependence of the dHvA amplitude. The compari- 
son of dHvA frequencies and m* between the experiment 
and calculation is given in Table [Q The effective mass 
in the experiment is enhanced roughly twice compared 
with the calculated one, which is consistent with the en- 
hancement of the Sommerfeld coefficient (7) between the 
experiment and the theory^ Such a large mass enhance- 
ment in the La compound is also found in LaFe4Pi2>^ 
The large mass enhancement in LaRu4Pi2 should be as- 
cribed to the electron-phonon interaction, taking into ac- 
count the rather high value of the superconducting tran- 
sition temperature Tq- Following the McMillan formula 
with the value of Coulomb pseudopotential ^i*=0.I3, 24 
the Debye temperature B D =405 K, and T c =7.3 
the electron-phonon coupling constant is estimated to be 




FIG. 5: The Fermi surface of LaRu4Pi2. The estimated vol- 
umes of the r- and 77-band FSs are ~ 0.3% and ~ 0.2% of BZ, 
respectively. 



0.65 in LaRu4Pi2, which is considered as a moderate- 
coupling superconductor. The value of 7 =44.4 mJ/mol 
K 2 and the jump of specific heat at Tc yield the ra- 
tio AC/jTq — 1.58, which is larger than 1.43 of weak- 
coupling BCS theory. 

The most probable reason why some of the theoreti- 
cal branches are not observed in the present experiment, 
particularly for the large dHvA orbits, is a strong re- 
duction of the dHvA signals due to their large curvature 
factors A" (A" = \d 2 A/dk^\, where A is the extremal 
cross-sectional area of the FS and fcn is the wave-vector 
component along the field direction)^ The rapid change 
of A as a function of fcn, i.e., the large value of A", dimin- 
ishes the dHvA amplitude for that extremal area. For the 
field angle 9 =30 degrees: A" = 0.14 for the a branch 
whereas A" = 22 for the ip branch in the calculation 
[see Fig. EJb)]. However, that is not the case for r and 
77 branches, which originate from 47 th - and 46th-band 
FSs, i.e., for H||[lll], A" = 0.18 for the r branch and 
A" = 0.61 for the 77 branch. The experimental limita- 
tions, i.e., insufficient temperatures and magnetic fields, 
could not be the origin of no observation of the r and the 
77 branches, since the [i branch with almost the same cur- 
vature factor (A" = 0.4) and the larger effective mass (see 
table[IJ than those of the r and 77 branches (the calculated 
effective masses for if||[lll] are 0.34mo for the r branch 
and 0.23too for the 77 branch) has been clearly observed. 
Therefore, dHvA signals should be observed if the t and 77 
branches really exist. The disagreement between the ex- 
periment and calculation concerning the r and r\ branches 
may be ascribed to the small difference of the lattice pa- 
rameter and/or the effect beyond the LDA treatment. 
Note that the band structure calculation shows the ab- 
sence of FSs responsible for the r and 77 branches in 
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TABLE I: Comparison of the dHvA frequency F and cyclotron effective mass m* between the experiment and the band structure 
calculation in LaRu4Pi2. 

Experiment Theory 
Field direction Branch F(xl0 7 Oe) m*/mo T(xl0 7 Oe) m*/mo 
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PrRu4Pi2, though the difference of lattice parameter be- 
tween LaRu 4 Pi2 (8.0608 A) and PrRu 4 Pi 2 (8.0424 A) is 
quite small (~ 0.2%) »2i Thermal expansion measurement 
suggests that the lattice constant decreases about 1.2% 
at low temperatures (<20 K) in LaRu^P^^ The band- 
structure calculation using the low-temperature lattice 
parameters will be clarified in future work. 

From the present dHvA experiment and the band- 
structure calculation, we have confirmed that the shape 
of the main FS of LaRu4Pi2 is like a distorted cube [see 
Fig. |5£a)] and its volume is nearly one-half of the first 
BZ, indicating the presence of good nesting instability 
with q =(1,0,0). The band calculation in PrRu 4 Pi2 also 
predicts the nearly perfect FS nesting, leading to the M- 
I transition>iS*ii The obtained topology of the main FS 
in LaRu 4 Pi2 is similar to that by band calculation in 
PrRu 4 Pi 2 . 

Then such a M-I transition should be expected also 
in LaRu 4 Pi2, however none has been reported yet; at 
least the resistivity and specific heat measurements 
show no M-I transition above Tq£& One can infer 
that the presence of two small spherical r and rj band 
FSs predicted by the band calculation, which are 
absent in PrRu 4 Pi2, might suppress M-I transition in 
LaRu 4 Pi2. However, the absence of r and rj branches 
in the present experiment rules out such a scenario with 
large reliability as mentioned above. In PrRu 4 P 4 2, the 
change of crystalline structureii*i2ii£ and the marked 
increase of resistivity below Tmi (Ref. 2) are explained 
by the nesting model, suggesting the disappearance of 



the entire FS below Tmi, based on the band-structure 
calculations with the localized 4/ electrons Note 
that in PrFe 4 P 4 2, which shows an antiferro-quadrupolar 
transition below 6.5 K, where only the main part of the 
FS disappears, the metallic state prevails and heavy 
fermion behavior appears&i^ The inelastic neutron 
scattering (INS) experiment in PrRu 4 P 4 2 shows sharp 
crystal electric field (CEF) excitation peaks below 
Tmi that gradually broaden across and above Tviir^ 
suggesting that the 4/ electrons are basically localized, 
though their hybridization with the conduction electrons 
(c — / hybridization) increases to some extent across 
and above Tvii- Cosidering all these facts, naturally, 
the present experiment suggests an essential 4/ electron 
contribution, which is the main difference between the 
two compounds, for the M-I transition. The c — / 
hybridization might cooperate with the FS-nesting 
condition in causing the M-I transition in PrRu 4 Pi2, 
which explains the absence of the M-I transition in 
LaRu 4 Pi2 with no 4/ electron contribution. 
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